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ABSTRACT 

^ \ By evolution of fermion mass matrices of the Fritzsch and the Georgi- 

an 

Jarlskog forms from the supersymmetric grand unified scale, DHR obtained 
predictions for the quark masses and mixings. Using Monte Carlo methods we 
test these predictions against the latest determinations of the mixings, the CP- 
violating parameter e K and the Bj-B® mixing parameter r^. The acceptable 
solutions closely specify the quark masses and mixings, but lie at the edges of 
allowed regions at 90% confidence level. 



One of the outstanding problems in particle physics is that of explaining the fermion 
masses and mixings. In the Standard Model (SM) the 6 quark masses, 3 charged lepton 
masses, the 3 quark mixings and the CP-violating phase of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix are introduced as phenomeno logical parameters. Over the years various mod- 
els have been proposed to reduce the number of these free parameters (JT) , of which the best 
known is the Fritzsch model f2|. Recently Dimopoulos, Hall and Raby (DHR) have pro- 
posed an ansatz for fermion mass matrices H] in the framework of minimal supersymmetric 
(SUSY) Grand Unified Theories (GUT). The DHR approach is based on the observation 
that some discrete symmetries present at the grand unification scale are broken in the low- 
energy theory. Thus some elements of the fermion mass matrices that vanish at the GUT 
scale are non-zero at the electroweak scale, and their low-energy values are calculable from 
the renormalization group equations. The fermion masses and mixings at the electroweak 
scale can thereby be expressed in terms of a smaller number of input parameters at the 
GUT scale. DHR work in the massless neutrino limit and relate the 13 SM parameters and 
a SUSY parameter tan/3 (discussed below) to 8 input parameters, leading to 6 predictions 
that include an allowed range of 147-187 GeV for the top quark mass (m t ). In comparison 
the Fritzsch approach gives 77 < m t < 96 GeV [|TJ, which is nearly excluded in the SM by 
the CDF experiment @ at a 90% confidence level (C.L.). 



The DHR quark mass matrices at the scale m t are 
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where all the parameters are real, tan/? = V2/V1 in terms of the Higgs doublet vacuum 
expectation values, and v = 246 GeV. The charged lepton mass matrix A4 e is obtained from 
the above form of M d by the substitutions = 0, 5 d = 0, E — > -3E', D —> D', F —> F' . 
At the SUSY-GUT scale, the parameters 5 U and 5d vanish and D = D', E = E' , F = F', 
so the input mass matrix A4 U is of the Fritzsch form and Aid and M. e are of the Georgi- 
Jarlskog form [f|, giving the GUT scale mass relations m;, = m T , m s ~ m^/3, ~ 3m e 



between quarks and leptons. The mass ratio prediction 



{m d /m s ){\ - m d /m s ) 2 = 9(m e /mJ(l - ra e /raj" 



(2) 



holds at all scales. 

The Wolfenstein parameterization || of the CKM matrix determined from the unitary 
matrices that diagonalize the DHR mass matrices can be expressed in terms of four angles 
(Oi) and a complex phase (0) as follows 



A = (sj + s\ + 2siS 2 COS 0) a = \V cd \ 



IK, 



X 2 A = s 3 — s 4 



iKfcl , 



AVp 2 + V 2 = s 2 = \V ub /V cb \ = ^m u /m c 
1] = Sis 2 sin0/A 2 , 



(3a) 
(3b) 
(3c) 
(3d) 



with Si = sin#j, q = cos^ (i = 1,2,3,4), where 6 2 , O3 are the angles that diagonalize the 
matrix Ai u , and 6\, 64 are those for Ai d @; only three of these angles are independent. 
These mixing angles are related to the quark masses and other parameters by 



Si ~ \jm d /m s , s 2 ~^Jm~Jrn c , s 3 ~\B/A\, s 4 ~ s 3 — |V cfc | . 



(4) 



The evolution based on the SUSY renormalization group equations (RGE) from the GUT 
scale to the appropriate fermion mass scales, taking all SUSY particles and the second Higgs 
doublet degenerate at the scale of m t M, gives the following relations, 
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x = (a G /a 1 ) 1/6 (a G /a 2 ) 3/2 , y = x{m b /m T )rj 1 / 2 rj b 1 



v(t*) = n ("gM 
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The RGE parameters 6j , q are given in Ref. 0. In these equations the couplings a± and a 2 
are evaluated at the scale m t . The mass parameters are defined as m q (fi = m q ) for quarks 
heavier than 1 GeV, and the lighter quark masses m s , ma, m u are calculated at the scale 
fi = 1 GeV. 

Starting from the well-determined values [0, a x {M z ) = 0.016887, a 2 (M z ) = 0.03322, 
and evolving at one- loop level to their intersection determines the GUT scale M G — 1.1 x 10 16 
GeV and the GUT coupling constant a G = 1/25 A. Evolving backwards, the strong coupling 
constant a s {Mz) = 0.106 is obtained, consistent with the LEP result a s = 0.118 ±0.008 |§. 
Also the values ai(mt) = 0.017 and a2(m t ) = 0.033 are determined, as well as the factors 
rj(m t ) = 9.7 and I(m t ) = 110. We have used a top-quark threshold of 180 GeV in the RGE, 
consistent with our output determination. In evolution below the electroweak scale we 
include 3-loop QCD and 1-loop QED effects in the running masses to obtain the evolution 
factors rjb = 1.44, i] c = 1.80 and T] s = 1.95 where rj q = m q (m q )/m q (m t ) for q = b,c 
and r] s = [m a (l GeV)/m M (l GeV)] j [m s (m t ) /m M (m 4 )]. Quark and lepton thresholds were 
handled by demanding that the couplings and running masses be continuous. The number 
of active flavors in the /3-functions and in the anomalous dimensions was changed as each 
successive fermion was integrated out of the theory. 

Following DHR, we take the following 8 relatively better-known parameters as in- 
puts: m e , m^, m T , m c , rrib, m u /md, |Kb| an d \V c d\- We generate random values for 
all inputs within 90% C.L. (1.64cr) ranges. The input mass values [f|,|nj are rn T = 
1784.1 + |^MeV, m c (m c ) = 1.27 ± 0.05 GeV, m b {m b ) = 4.25±0.1GeV, where la er- 
rors are quoted. We also impose the theoretical constraint 0.2 < m u /md < 0.7 JTT). We 
next calculate rrid, m s from Eqs. (0) and (|5aD (obtaining = 5.93 MeV, m s = 146.5 MeV), 
Sx and s 2 from Eq. (^), s 3 from Eq. ( pb]) for the input of \V c b\, S4 from Eq. (^), and (ft from 
\Vcd\ of Eq. ([3a]). Using these values we evaluate the magnitudes of all elements of the CKM 
matrix. We retain only those Monte Carlo events that satisfy the following ranges from the 
1992 Review of Particle Properties [EJ, 
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CKMl 



0.9747-0.9759 0.218-0.224 0.002-0.007 
0.218-0.224 0.9735-0.9751 0.032-0.054 
0.003-0.018 0.030-0.054 0.9985-0.9995 



(7) 



as well as the ratio 0.051 < \V u b/V c b\ < 0.149 ||. Figure [I] shows a scatter plot of sin/3 
versus m t obtained from our Monte Carlo analysis; one sees that only a narrow wedge of 
the space is permissible. 

From sin/3 < 1 in Eq. (5), \V c b\ must satisfy the inequality 

.1/2 
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31 



m b m c 



> 0.053 



(8) 



ttv V 2(1 — y 12 ) ■m T r\r]\ ) r] { 

which is just at the edge of the 90% C.L. allowed range. Calculating m t from Eq. (|5a|) and 
sin/3 from Eq. (|5H) and requiring that \V c b\ be within its allowed range, we find 



174 <m t < 183 GeV , sin (3 > 0.954 (tan (3 > 3.2) 



(9) 



This top quark mass determination is consistent with estimates from the electroweak radia- 
tive corrections || but is much more restrictive. The predicted value of tan j3 is large, which 



may have significant phenomenological implications for Higgs boson searches at colliders |12 
Next we include the constraints from the measured values 



~K\ 



(2.259 ± 0.018) x 10" 3 [9] , r d = 0.181 ± 0.043 @ , 



(10) 



of the CP- violating parameter e K and the B^-B^ mixing parameter r^, The theoretical 



formulas, including QCD corrections, can be found in Eqs. (2.1) and (2.10) of Ref. [p~4 1 . In 
our Monte Carlo analysis we allow variations of the bag-factors and -B-decay constant over 
the following ranges fT4|j : 



0.33 < B K < 1.5 , 0.1 GeV < JB B f B < 0.2 GeV 



(11) 



taking f K = 160 MeV and AM K = 3.521 x 10~ 15 GeV. The solutions so obtained closely 
specify the CKM matrix to be 
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^ 0.9749-0.9759 0.2185-0.2230 0.0027-0.0032 ^ 



Wt 



CKMl 



0.2185-0.2230 0.9735-0.9745 0.0530-0.0540 
0.0106-0.0109 0.0518-0.0529 0.9985-0.9986 



(12) 



and predict the CP-violating phase to be in the range 



70° < 6< 80°. 



(13) 



The inclusion of e K and r d almost uniquely determines the values of \V t d\ and \V ts \. Since 
| Vcb\ is near its allowed upper limit, \V u b\ is pushed to its lower end by the unitarity condition. 
The output value of the ratio 



0.051 < \V ub /V cb \ < 0.059 



(14) 



is at the low end of the allowed range. Improved experimental determinations of \V c b\ and 
\V u b/V c b\ will test the DHR ansatz. In terms of the Wolfenstein parameters 0, we find 



0.2185 < A < 0.2230 , 1.07 < A < 1.13 , 
0.195 < T] < 0.243 , 0.105 < p < 0.129 , 0.222 < Vp 2 + V 2 < 0-275 

The output values of the mass ratios of the light quarks are 

0.52 < m u /m d < 0.70 , m d /m s = 0.0405 , 0.021 < m u /m s < 0.028 



(15) 



(16) 



giving 3.08 < m u < 4.15 MeV. These light quark masses and their ratios are consistent with 
those obtained in Refs. [0,01 > do no ^ a g ree as wen with some other recent studies |15 
in which m u /md 0.3 was obtained. 

The heavy quark masses are now constrained to the narrow ranges 



1.19 <m c < 1.23 , 4.09 < m b < 4.20 . 



Another interesting result is restrictive ranges for the constants Bk and Jb 



(17) 



0.33 <B K < 0.43 , 0.14 < J B B f B < 0.17 GeV 



(18) 
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on which theoretical uncertainties have been problematic |16 



We conclude with some brief remarks. From Eq. (pig) , m t is inversely proportional to 
VbVcV 1 ^ 2 \ an d the theoretical uncertainty in this quantity could somewhat enlarge or close 
the window in m t (and correspondingly the window in |V^,|). The DHR analysis assumes 
dominance of the top quark Yukawa couplings in the RGE evolution. Since the output tan /3 
may be large, the effects of fully including and A r in the evolution may not be negligible; 
this question deserves further study. Also two-loop renormalization group equations between 
Mz and Mgut should eventually be incorporated. We have studied the charged Higgs boson 
effects on e K and r^. With M#± degenerate with m t , as assumed in the model, we found no 
significant changes in our results. This is due to the fact that H ± effects are smaller at large 
tan p for the K and B systems. In summary, the DHR ansatz for fermion mass matrices is 
consistent with all current experimental constraints at 90% C.L.. It leads to almost unique 
values for m t and quark mixings which make it an interesting target for future experiments. 
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FIGURES 

FIG. I. Scatter plot of sin (3 versus m t from our Monte Carlo analysis of the DHR model, 
imposing the constraints of input masses and present values of CKM matrix elements. 
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